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Broad Range of Flow Problems of Interest 

Wide Range of Flow Conditions: 

Subsonic - Hypersonic 

Internal - External - Store Separation 

Cruise - High Angle of Attack 

Rows phenomena of Interest: 

Inlets/Diffusers 
Streamwise Curvature 
Shock/BL Interactions 
Rectangular Duct — Circular 

Nozzles 

Entrainment 

Round— Rectangular Duct 
High Speed Shear Layers 

External Aerodynamics 

Vortex 

Leading Edge Separation 
Shock/BL Interactions 


Leading Edge Separation - Cowl Ups 
Separation Induced Unstart 


Film cooling, Uners, Vanes 
Swirl 


3D Boundary Layers 
Wakes 


The CFD Environment at Lockheed Fort Worth Company 


Most codes developed or highly modified in house 
General grid generation and solvers for diverse applications 
Structured and unstructured solvers 
Computational efficiency important 

• Complex geometries, many gridpoints 

• Large arrays of flow conditions 
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Requirements for Turbulence Models 


Turbulence Modeling Priorities for Industrial Application 

• Validation 

High accuracy for attached flows 
Reasonable accuracy for all flows 
High confidence level 

• Computational efficiency 

• Robust for complex geometries 

• Transitional modeling capability 


To 


obtain acceptable accuracy, propulsion flows demand more sophisticated 
turbulence models than do external aerodynamic flows 


The k - kl and k - 1 Two Equation Turbulence Models 

Advantages of using kl or I Instead of e or to 

kl and I equations are easier to resolve numerically than e equation 

Dissipation Length Scale is an integral length scale 

•Can derive equation for volume integral of two point correlation function. 
•Theoretical e equation Is dominated by small scales 
k - kl and k - 1 agree better with compressible boundary layer data than 
does k-z 

Disadvantage - current formulation requires calculation of distance to walls 


k - kl model 

• Includes unique, consistent wall 

function 

• Accurate for transonic flows 


k - \ model 

• Derived from k - kl model - Identical in 

high Re turbulence 

• Near wall model simulates k In viscous 

sublayer 
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The k - kl Model Wall Function 


Wall layer model derived from and consistent with the k - kl model 

• Assume convection In momentum, energy and turbulent kinetic energy 

equations to be negligible 

• Boundary layer approximation 

Match velocity, k and I at first grid point In Navier - Stokes solution 
First grid point can be in viscous sublayer, buffer or logarithmic region 
Boundary conditions on k and I simple for k - kl model 
Advantages of wall functions 

• Reduces number of necessary grid points 

• Reduces number of iterations to converge steady state solution 60 - 90% 


Wall Functions are Accurate for Separated Flow Applications 




Axisymmetric Bump, Transonic Flo w Experi ment 


Accurate predictions with and without wail 
functions 




Velocity profiles with and without wall functions 
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The k - 1 Model with Near Wall Model 

kl equation is transformed exactly to an I equation 

Advantages of k - 1 formulation 

• I is linear near wall, kI nonlinear and very small 

• Near wall damping terms disappear 

• Production term drops out with current choice of constants 

k - 1 model includes: 

• Transitional flow modeling 

• Compressibility corrections 

Modeling of details of k profile near wail important for hypersonic flows 

• Magnitude of normal stress term comparable to static pressure 

• Near wall density variations large 


/ Equation Much Easier to Resolve than e Equation 

€ equation requires fine grid from wall to y + of 20 to resolve peak 

•Exclusion of near wail viscous dissipation term aggravates problem 
•Logarithmic region, e« 1 /y 

I equation Is nearly linear near wall - much less sensitive to grid resolution 
Length scale and dissipation profiles near wail 
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Resolution Study with k - e and k-l Models 
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k - e with Lam * Bremhorst 



Symbol 

Number ol 
Grid 
Points 

Stretching 

Rate 

from wall 

Y*, first 
grid point 

— 

151 

1.04 

0.033 

□ 

11 

1.4 

4.3 

0 

11 

1.6 

1.65 

A 

11 

1.3 

0.67 

0 

11 

1.9 

0.44 


Sample Applications: 

Mach 8 Shock Wave Turbulent Boundary Layer Interactions 
F-16 Inlet Derivitive, Isolated Duct Study 
Multi-slot Ejector 
F1 10 Nozzle Drag Reduction Study 
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k - I Model With Compressibility Correction gives Best Prediction 
For Mach 8 Shock Boundary Layer Interaction 





The k - I Model Predicts Turbulent Shock - Wave Boundary 
Layer interaction Well 



Mach 8, 10 Degree Wedge Generator 
2D case, Separated Row 



Fin* Grid Solution, 107x101 
— Com* Grid, 97x111 
O A Experiments* Oat* 
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Good Predictions of Multi - Slot Ejector Obtained with 

k - kl Model 
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Mach Contours 


k - kl Model Predicts Entrainment Effects Near Slots 


Velocity vectors colored by Mach Number 
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Summary 


Computationally efficient k - 1 and k - kl models have been developed 
and Implemented at Lockheed Fort Worth Company 


Many years of experience applying two equation turbulence models to 
complex 3D flows for design and analysis 


37 




